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Abstract 



^ c*| We show that the more energetic superluminal neutrinos with quadratically dispersed 
P-l superluminalities 5 = (3 2 — 1, for /3 — v/c where v is the neutrino velocity, also lose 
P_| significant energy to radiation to the v + e~ + e + final state in traveling from CERN to 
Gran Sasso as has been shown to occur for those with constant superluminality by Cohen 
and Glashow if indeed 5 ~ 5 x 10~ 5 . In addition, we clarify the dependence of such 
radiative processes on the size of the superluminality. 



In trying to understand the results of the OPERA Collaboration [T] for apparent 
faster than light velocity for neutrinos in transit from CERN to Gran Sasso, Cohen and 
Glashow [2] have made an estimate of the Cherenkov effect to the v + e + + e~ final state 
using the standard methods adapted to Lorentz invariance- violating interacting quantum 
field theory. This effect is also studied in Ref. [3] in a general framework for such violations 
of Lorentz invariance with particular attention to the case where the superluminality is 
quadratically dispersed in contrast to the case considered by Cohen and Glashow with 
constant superluminality, where we define the superluminality by 5 = (3 2 — I when v is the 
neutrino velocity. Here, we extend the analysis of Cohen and Glashow to the quadratically 
dispersed case of Ref. [3] and we compare the analyses of Cohen and Glashow and of 
Ref. [3] with an eye as well toward the very interesting observations of the OPERA 
Collaboration. 

We need to first state very clearly at the outset here that we do not claim to rule out 
the possibilities discussed in Refs. |IH6] in which appropriate deformation of the Lorentz 
symmetry renders these Cherenkov processes kinematically forbidden just as one has for 
the usual non-superluminal case in the usual special theory of relativity. Thus, non- 
observation of the effects discussed herein does not rule-out such changes in the standard 
orthodoxy. 

Specifically, Cohen and Glashow arrive at the results 

r = k '-^-E 5 5 3 

1927T 3 /-Q 

dx 192tt 3 

for the Cherenkov radiation rate and energy loss per unit distance for the process v — > u+ 
e + +e~ for u's of energy E where the parameters k', k are 1/14 and 25/448, respectively [7J. 
In this way, assuming that 6 is independent of E, they find [2] that the z/s should lose 
considerable energy between CERN and Gran Sasso. 

On the other hand, Mattingly et al. in Ref. [3] have discussed the same transition 
for superluminal i/'s via analogy with the so-called //-splitting process v\ — > v a v b^b for 
which they find the result 

VVD ~ 64vr 3 M Pl [ ] 

where the factor of 3 multiplying Gp sums over three neutrino species so, ignoring the 
masses of the neutrinos and the electron, Mattingly et al. assert (see Sect. 6 of Ref. [3]) 
that the rate for the v + e + e~ decay should be nearly the same as that for the ^-splitting 
decay in fl2}T_!]. The Lorentz-violating parameter is such that the dispersion relation for 
the v is 

E 2 p =ml+p 2 + r ]u p' l /M 2 Pl 

where Mpi is the Planck mass, m„ is the neutrino rest mass and p is the v momentum. 
The parameter r\ v is to be determined from experiment. Thus, using v = and ignoring 
the neutrino mass we see that 

5 ^ 3r] u E 2 /M Pl 



1 The sum over the three neutrino species mildly over compensates for the almost equal coupling of the 
Z° boson to left and right handed electrons so that the rates would become almost equal if we multiply 
latter one by 2/3 - we ignore such factors here. 
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so that Mattingly et al. argue that the rate for the v — > v + e + + e process should be 
nearly given by 



r + - ~ r - 

1 ue + e — 1 vvv 
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64vr 3 Af|, (3) 



64(81)tt s 

While the numerical factors and the explicit dependence on E are consistent between the 
two estimates of T for the v — > v + e + + e~ process, the two dependencies on 5 differ by 
one factor 5: this factor would obviate the arguments of Cohen and Glashow, given that 
OPERA have found that 5 = 5 x 10 -5 . Which dependence on 5 is actually correct? The 
literature is beginning to feature use of both dependencies, as we see in Ref. [6] where 
the corresponding dependence of Mattingly et al. is used and in Ref. [8] where both 

dependencies are featured Q. Secondly, what effect does quadratic dispersion for 5 have on 
the results of Cohen and Glashow? We address these two questions here in inverse order. 

We see that, unlike the case discussed by Cohen and Glashow, we have to consider 
pronounced dependence of 5 on momentum: for a neutrino of initial energy Eq ~ 20GeV, 
5(E) = 5q(E/Eo) 2 in an obvious notation. Moreover, the factor of 3 between the devi- 
ation between the squared group velocity and 1 here and the squared incoming v four- 
momentum vector in units of its squared laboratory energy E means that we should 
identify the quantity 5 in the Cohen-Glashow formulas in ([T]) with our 5/3. Thus, the 
dependence of dE/dx on E here is given by 

jg ., G ? E " S3 » (4) 
dx 192tt 3 ES ' W 



Integrating this relation gives us 
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where we introduced the distance L = 730km from CERN to Gran Sasso and the Cohen- 
Glashow quantity \ jE\ = 5A g|j L so that for 5 = 5 x 10" 5 E T = 12.8GeV. This should 
be compared with Eq.(21) in Ref. [8], wherein only one power of the analog of i] v appears 



2 By going to the 'preferred frame' of the incident v, we see from the textbook results in Ref. [S] 
that, when we neglect all final state particle masses, the total rate to the ve~e + final state is in fact 
independent of the handedness of the outgoing electron when v e , the SM vector coupling of e~ to the 
Z°, is set to zero. Thus the arguments in Ref. [8] as to why their splitting process rate has a different 
dependence on S arc not valid. 
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versus three powers here. If we define the average value of 5 via0 



< 5 > = — / dxS(x) 
L Jo 




[ dx(E 2 (x)/E 2 )S , 



(6) 



we arrive at the relation 
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This means that, if we identity 3 < 5 > with the OPERA result 5 x 10~ 5 , we can solve 
this latter relation (J7|) to get 3<5o = 5.3 x 10~ 5 when the initial energy is E = 20GeV. 
Thus, even though 5 varies quadratically with momentum, it does not deviate very much 
from its initial value and a 20 GeV v now loses about 5.8% of its energy in travelling from 
CERN to Gran Sasso compared to 37.4% for the case discussed by Cohen and Glashow. 
The two groups of average velocity data in Ref. [1] have average energies 13.8 and 40.7 
GeV with superluminalities of - 1 given by (2.25 ± .81) x 10" 5 and (2.791 jg) x 10~ 5 
respectively. Solving for the respective values of 3<5o we get 5.1 x 10~ 5 and 1.2 x 10~ 4 with 
energy losses of 1.1% and 41.8% respectively, where we use again 3 < 5 >= 5 x 10~ 5 in 
view of the errors on the data without loss of content - the respective energy losses in 
the Cohen-Glashow analysis are 16.6% and 68.6%. The factor of order 2 variation of Sq 
with initial creation energy is not excluded by the data. In any case, when the v energy 
is in the 10 MeV regime, our values of 5 are well below the limits from SN1987a [TO] . 
4 x 10~ 9 , due to the quadratic dependence on energy in principle. We say 'in principle' 
because there is, presuming the OPERA result is true, an unknown source of this energy 
dependence and possibly of the variation of 5q with initial energy and we cannot make 
definitive statements until this source is known. We need more experimental input. 

Concerning the detailed relationship between our calculation and that in Ref. [TTJ , we 
note that the authors in Ref. [TTJ make a numerical estimate for the v — > v + e + + e~ 
process with the quadratic dispersion for the v superluminality and with the usual special 
relativity dispersion for the e~, e + with the result that the coefficient of G 2 F E 5 5 3 (1 — 
4sin 2 6 W + 8 sin 4 6 W ) in the Cohen-Glashow process is 1/45197.8 vs the exact result [2JE] 
1/41672.4, neglecting the lepton rest masses. This is acceptable in view of the current 
errors on the data. The authors in Ref. [JTJ then approximate the energy loss using a 
constant value of 5 so that they give in their eq.(21) a mean decay time that scales as 
the 5th power of the v energy when in fact the dependence of 5 on energy shows that 
the decay time has to scale with the 11th power of this energy, as we have shown here. 
This results in significant deviation from the correct decay energy loss profile as we have 
shown above for the higher energies in the OPERA spectrum. However, for the 17.5 GeV 
energy z/'s, we find that they lose 3.2% of their energy in travelling from CERN to Gran 
Sasso whereas the authors in Ref. [TTJ estimate this as a 4.2% loss if you take their result 
literally as written below their eq.(20). If you interpret their 4.2% as the mean number 

that decay via the v — » v + e~ + e + process so that we do not actually lose z/'s but we 
lose their energy, then to get an estimate of the corresponding energy loss we have to 

3 Note that the strict definition of the average velocity, < v >— J dtv(t)/T = Lj J dx/v(t), reduces to 
this result in ([6]) in the limit in which we work wherein all velocities are close to 1. 
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multiply by the mean energy loss per decay, which is j2j[7] 78%. This would then agree 
with our 3.2% result - this is expected because as the energy does not change very much 
the assumption of a constant 5 is not a bad one. We note as well that the authors in 
Ref. [IT] do not address the issue of the difference between the power of 5 in the decay 
rate formulas in ([1]) on the one hand and in (T5]) and ([3]) on the other. We now turn to 
this issue. 

To answer the question on the proper power law for the dependence of the radiation 
rate on 5, we revisit the analysis in Ref. [3] and observe that, in their Eq.(4.10), there 
is an energy conserving delta function, 5(E P — E' p — E q — E' q ), where for definiteness we 
reproduce this equation here as 



r 



AB 



9' 



j 5(E P -E' p -E q - E' q )d 3 p'd 3 qF 2 



16(2tt) 5 M|cos 4 ^ 

for the width for the va(p) 1/ a{v') 1/ b{. ( i')i / b{. ( 1) splitting process in a standard notation, 
with E' a = E a >, a = p, q. Here, q'(q) is the 4-momentum of the z/ B (z/ s ), respectively, 
and the SU 2 l xUi [12] couplings and parameters are standard by now - the rest mass of 
the Z° is M z , the SU 2 l coupling constant is g, the Weinberg angle is 6 W . The function 
F contains the dependence of the differential decay rate on the various angular phase 
space variables [3] . Upon arguing correctly that the excess energy from Lorentz- violation 
goes to transverse motion, Mattingly et al. then estimate the size of this transverse 
phase space and argue that the only effect of the delta function constraint on energy is 
to remove one factor of the energy from their result for the integral over the two phase 
spaces for p' and q. However, there is a Jacobian required to implement this removal, by 
the standard methods, given by the absolute value of the derivative of the argument of 
the delta function with respect to pj, for example^ due to the result 

*(/(*))= E ( 9 ) 



x a (z{zeras of /} 



I dfjXg) 

dx 



for any sufficiently well-behaved function f(x) : This provides the missing factor of S that 
should divide the expressions on the RHS of ([3]) and of ([2]) to bring them into agreement 
with Cohen and Glashow. To see this note that the absolute value of the derivative of 
the argument of the delta function with respect to pj, can be represented as follows: 
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6(1 — x — y)< 



Here, || denotes the direction along that of the 3-momentum of the initial state v. 
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where we follow Ref. [3] and use the longitudinal momentum fractions x, y for p' and q 

2 

respectively and where we parametrize the transverse momenta as r_j_ = f±y/rj^^—\ r for 

the unit vector f± in their direction. By going to the "preferred frame" of the incoming v 
we see that the parameters 2A a , defined to be non-negative here, are bounded by constants 

close to 1. We have used the usual dispersion relations E q = \/ q 2 + mf, E' = \/ q' 2 + m\ 
for the e + and e~ respectively, where m e is the rest mass of the electron and have neglected 
all electron and neutrino massed. We stress that this factor of 5 in ([TO]) should divide 
the RHS of ((21) also when one uses it for the process v —> vvv by the standard methods. 

We conclude that the estimates of Cohen and Glashow that the very energetic z/s 
should lose a large amount of their energy due to Cherenkov e + e~ radiation hold true for 
the observations of OPERA in Ref. [I] for superluminal u's with quadratically dispersed 
values of S as well. While showing this result, we have clarified the issue of the dependence 
on 5 of this and the attendant related radiative processes. 

Note Added in proof. - We understand [13] that the difference noted in the text 
between our Eq.(J3]) and Eq.(21) in Ref. [8] did not enter the respective results given in 
Ref. [8]. 
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